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■ ABSTRACT 

' We present a method to recover mass profiles of galaxy clusters by combining data on thermal 

, Sunyaev-Zeldovich (tSZ) and X-ray imaging, thereby avoiding to use any information on X- 

ray spectroscopy. This method, w hich repr esents a development of the geometrical deprojec- 
tion technique presented in Amegho et alj (i2007). implements the solution of the hydrostatic 
^ . equilibrium equation. In order to quantify the efficiency of our mass reconstructions, we ap- 

Q>^ ' ply our technique to a set of hydrodynamical simulations of galaxy clusters. We propose two 

Q>^ I versions of our method of mass reconstruction. Method 1 is completely model-independent 

and assumes as fitting parameters the values of gas density and total mass wit hin different 
CsJ ' radial bins. Method 2 assumes instead the analytic mass profile proposed by .Navarro et al.l 

(Il997h (NFW). We find that the main source of bias in recovering the mass profiles is due to 
deviations from hydrostatic equilibrium, which cause an underestimate of the mass of about 
' 10 per cent at r^oo and up to 20 per cent at the virial radius. Method 1 provides a reconstructed 

. mass which is biased low by about 10 per cent, with a 20 per cent scatter, with respect to the 

true mass profiles. Method 2 proves to be more stable, reducing the scatter to 10 per cent, but 
. ^ ^ with a larger bias of 20 per cent, mainly induced by the deviations from equilibrium in the 

. outskirts. To better understand the results of Method 2, we check how well it allows to recover 

' the relation between mass and concentration parameter. When analyzing the 3D mass profiles 

5^ I we find that including in the fit the inner 5 per cent of the virial radius biases high the halo 

concentration, thus suggesting that the NFW profile is not a perfect fit in the central regions 
of our simulations including cooling and star formation. Also, at a fixed mass, hotter clusters 
tend to have larger concentration. Our procedure recovers the concentration parameter essen- 
tially unbiased but with a scatter of about 50 per cent. In general, our analysis demonstrates 
that combining X-ray imaging with spatially resolved tSZ data is a valid alternative to using 
X-ray spectroscopy to recover the mass of galaxy clusters. 

Key words: large-scale structure of Universe - galaxies: clusters: general - cosmology: mis- 
cellaneous - methods: numerical 



1 INTRODUCTION 



A number of important cosmological tests are based on mass mea- 
surements in galaxy clusters. In particular the mass function, the 
baryon fraction and their redshift evolution are highly sensitive to 
the underlying cosmology and provide constraints on the Dark Mat- 
ter and Dark Energy content of the Universe. Precise mass mea- 
surement in galaxy cluster are then necessary to calibrate clusters 



as precisio n tools for cosmology ("Haiman et al."200l'; iRosati et al.l 
[2002 ; Pierpaoli et al.ir2003. : .Voit.2 005: Borgani 2006). 

In X-ray studies, the total collapsed mass of a cluster is 
determined by applying the hydrostatic equilibrium equation to 
gas density and t emperature prof i les. A number of authors (e.g., 
iRasia et alj|2004l ; iKav et aLll2004 iNagai et alj|2007l) analyzed hy- 
drodynamical simulations of galaxy clusters and found that the 
gas is not in perfect hydrostatic equilibrium. Instead, they found 
deviations up to 20 per cent, which have the effect of systemat- 
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ically biasing low the reconstruction of the total collapsed mass. 
The amount of this underestimate depends on both the model 
assumed and the dynamical state of the cluster. For instance, 
iRasia et alj jZOOd) showed that ass uming the isothermal /3-model 



I Cavaliere & Fusco-Femianol 1976h fo r the g as distribution gives 
the worst reconstruction. iNagai et al] feOOTD differentiated their 
sample of simulated clusters in relaxed and unrelaxed objects, with 
the latter showing a large r scatter in the mass reconstruction. Fol- 
lowing the same direction, IJeltemaetai] (12007.) found a correlation 
between quantitative measures of the morphology of the X-ray im- 
ages and the bias in the mass reconstruction, although with a quite 
large scatter. lPuchwein & Bartelmannl ( l2007h probed the deviations 
from hydro static equilibrium in di fferent stages of a merger pro- 
cess, while IPiffaretti & Valdaminil fcoOS) performed an extended 
analysis of a large sample (~ 100) of simulated galaxy clusters, in 
order to disentangle various biases in the mass reconstruction. 

Recovering mass profiles via the hydrostatic equilibrium 
equation involves the derivatives of both gas density and temper- 
ature profiles. Then, an accurate mass determination requires high- 
quality temperature measurements. For this reason. X-ray studies 
are often limited to the inner regions of the clusters and to objects 
at moderate redshift. Furthermore, all the temperature profiles so 
far used for the reconstruction of cluster masses are bas ed on X- 
ray data. As ori ginally noticed by iMazzotta et al.l j2004) (see also 
IVikhlinirj|200^ . the thermal complexity of the ICM causes the 
temperature determined by fitting the X-ray spectrum to a single- 
temperature bremsstrahlung model to be generally different from 
the electron temperature. On the other hand, it is the electron tem- 
perature that enters in the equation of hydrostatic equilibrium (un- 
der the assumptions of fully ionized plasma and equilibration be- 
tween ions and electrons). The question then arises as to whether 
the difference between spectroscopic and electron temperature may 
induce an additional bias in the mass estimate. 

In this paper, we propose to use a combination of X-ray imag- 
ing and tSZ data, which avoids the use of X-ray spectroscopy. 
Taking advantage of the different dependence of the tSZ signal 
and of the X-ray emissivity on gas density and temperature, one 
can recover both by suitably deprojecting tSZ and X-ray imag- 
ing data, without the need of any X-ray spectroscopic informa- 
tion. Furthermore, both X-ray and SZ cluster images can be ob- 
tained in principle out to larger radii than possible for X-ray 
spectra. Therefore, their combination could allow studies of the 
ICM thermal properties to be pushed out to larger fractions of 
the whole cluster virialized regions. Indeed, different algorithms, 
which use a combination of tSZ and X-ray data, have been pro- 
posed by a number of authors and applied to analytical cluster 
models and/or sets of simulated clusters (e.g. 'Zaroubi et al."200l'; 
Lee & Suto 2004; De Filippis et al. 2005; Puchwein & Bartelmann 
20071 ; lAmeglio et al.ll2007 . and references therein). Owing to the 
relatively low angular resolution of currently available tSZ tele- 
scopes, this metho d has been so far applied only to a handful of 
observed clusters jZhang & Wiil 120001 ; iPointecouteau et a i] |2002l ; 



iMroczkowski et al.l2008l) . with results that are generally consistent 
with those based on X-ray spectroscopic data. 

Accurate mass profiles reconstruction can also be used as 
probes for cosmology. In fact, mass profiles are expected to fol- 
low a universal functional form, which is valid over a wide range 
of masses, from dwarf galaxies to massive galaxy clusters. A for- 
mulation for this function has been originally prov ided by fitting 
the ma ss profiles in a set of N-body simulations by iNavarro et al.l 
( 1 19971) (NFW hereafter). The analysis of X-ray data from Chandra 
and XMM-Newton observations confirm the validity of the NFW 



model, out to a significant fraction of the cluster v i rial radius (e.g. 
I Pratt & Aniaudl 12 002; Poi ntecouteau et al.l lioosi ; Ivikhlinin et al.l 
boOd lZappacosta et all200^ An important parameter of the NFW 
model is the concentration c, which is given by the ratio between 
the halo's virial radius and the characteristic radius of the density 
profile. In the hierarchical structure formation scenario the more 
massive objects are expected to form recently, from a lower-density 
environment than less massive ones. Thus, it is expected an in- 
verse correlation between the mass of an object and its concentra- 
tion, with a substantial scatter related to th e distribution of halo 
formation epochs (e.g., iNavarro et al.l [19971 ; iBullock et alj|200ll ; 
iDolag et al.ll2004h . This relation has been confirmed by the obser- 
vation of m ass profile s in gala xy c l usters, using both X-ra y and 
lensing data l lSchmidt & Allen 20071 ; iMandelbaum et al.|[20ba. and 
references therein). These authors generally agree in finding a well 
defined relation with a substantial intrinsic scatter, as expected from 
the predictions of numerical simulations. However, some discrep- 
ancies are still present in the determination of the slope and nor- 
malization of the relation, which may be generated by biases in the 
mass measurements and/or by selection effects. 

We ex tend here the deprojection algorithm presented by 
lAmeglio"et al. (2007) (A07 hereafter), by including the solution of 
the hydrostatic equilibrium equation. Using this technique, we an- 
alyze a set of 14 simulated clusters having Tsi^ 3 keV, with the 
aim of quantifying the accuracy with which total mass profiles can 
be recovered by combining X-ray and tSZ images. We will also 
discuss how the relation between halo concentration c and mass 
M can be recovered with our deprojection method, als o compar- 
ing it with the theoretical predictions of the model by lEke et al] 
( I2OOII) . The analysis presented in the following can be applied to 
data from the present generation of X-ray telescopes. Furthermore, 
future X-ray observations with lower background (e.g. as expected 
from the eROSITA missiorQ) will provide good imaging data for a 
large number of clusters out to 2; ~ 1. As for the tSZ data, exploit- 
ing the full potentiality of our technique would require spatially 
resolved observations, which will be available from upcoming (or 
just started) SZ experiments, based both on interferometric arrays 
(ALMA: Atacama Large Millimeter Arra50; CARMA: Combined 
Array for Research in Millimeter-wave Astronomjlfl) and on single 
dishes with large bolometer arrays (CCAT: Comell-Caltech Ata- 
cama Telescope|3; LMT: Large Millimeter Telescop43). 

The paper is structured as follows. In Section |2] we briefly 
present the set of simulated clusters. We use the same subset of 
simulated clusters having Tsif^ 3 keV, which is described by A07. 
In Section[3]we introduce the hydrostatic equilibrium equation and, 
as a preliminary test, probe the intrinsic deviations from this equi- 
librium for our set of simulated clusters. After briefly reviewing 
the deprojection algorithm introduced by A07, Section|4]describes 
the implementation of the hydrostatic equilibrium equation in this 
algorithm. Sections |5] and |6] present our results on the total mass 
reconstruction and on the c-M relation respectively. Finally, our 
main conclusions are summarized in Section|7] 
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2 THE SET OF SIMULATED CLUSTERS 

The sample of simulated galaxy clusters used in this paper has been 
extracted from a large-scale cosmological hydro-N-body simula- 
tion of a "concordance" ACDM model with f2,„ = 0.3 for the 
matter density parameter, f^A = 0.7 for the cosmological constant, 
fib ~ 0.019^"^ for the baryon density parameter, h = 0.7 for 
the Hubble constant in units of 100 km s^^Mpc^^ and as ~ 0.8 
for the r.m.s. density perturbation wi thin a top-hat spher e having 
comoving radius of 8/i~^Mpc (see 'Borgani et al. l l2004l . for fur- 
ther details) . The simulatio n, performed with the Tree-l-SPH code 
GADGET-2 jSpringell [20o3) . follows the evolution of 480^ dark 
matter particles and an initially equal number of gas particles in 
a periodic cube of size 192/i~^ Mpc. The mass of the gas particles 
is mgas = 6.9 X lO*/i~^M0, and the Plummer-equivalent force 
softening is 7.5h~^ kpc at z = 0. The simulation includes the 
treatment of radiative cooling, a uniform time-dependent UV back- 
ground, a sub-resolution model for star formation and e nergy feed- 
back from galactic winds jSpringel & Hernquislll2003h . At z = 
we extract a set of 117 clusters, whose mass, as computed from a 
friends-of-friends algorithm with linking length 6 = 0.15 (in units 
of the mean interparticle distance) is larger than 10^* h^^MQ. For 
these clusters we compute the spectroscopic-like temperature 



,a-l/2 



(1) 



where a = 0.75 is a parameter shown bv lMazzotta et alj ( |2004) to 
accurately reproduce the temperatures obtained from the spectro- 
scopic fit. 

Due to the limited box size, the largest cluster found in the 
cosmological simulation has T^i = 4.6 keV. In order to ex- 
tend our analysis to more massive and hotter systems, which are 
mostly relevant for current tSZ observations, we include four more 
galaxy clusters having i\fvir > 10^^ /i^^A/glfl and belonging to 
a different set of hydro-N-body simulations (Borgani et al. 2006i). 
These objects have been select ed from a DM-only simulation of 
a large cosmological volume ( lYoshidaetaIll200lh . and resimu- 
lated at higher resolution. The achieved resolution corresponds to 
JTigas = 1.69 X IQ'^ Mq for the mass of the gas particle and 
a gravitational softening of 5/i~^kpc at z = 0. These simulations 
have been performed with the same choice of the parameters defin- 
ing star-formation and feedback. The cosmological parameters also 
are the same, except for a higher power spectrum normalization, 
(78 = 0.9. 

From this large set of simulated clusters, we select a subsam- 
ple of objects having Tsi^ 3 keV, which are more relevant for SZ 
studies. We identify 15 suitable objects, from which we exclude 
one very irregular cluster. The main characteristics of this subset of 
simulated clusters are listed in Table [T] where the cluster labeled 
C4, C12, C13 and C14 are those based on the high-resolution res- 
imulations. All the results shown in this paper are based on the 
analysis of these 14 objects. A larger subsample of systems with 
Tsi < 3 keV, extracted from the same box, are used for the analy- 
sis presented in Section 6. 



° Here and in the following, the virial radius, rvir, is defined as the ra- 
dius of a sphere centred on the local minimum of the potential, containing 
an average density, Pvir> equal to that predicted by the spherical collapse 
model. For the cosmology assumed in our simulations at low redshift it is 
Pvir — lOOpc, being pc the cosmic critical density. Accordingly, the virial 
mass, Afvir, is defined as the total mass contained within this sphere. 
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2.0 


1.0 
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3.0 


2.0 


1.0 


Cll 


6.3 


3.0 


2.2 


1.0 


C12 


32.0 


8.9 


3.8 


1.9 


C13 


19.2 


6.3 


3.2 


1.6 


C14 


19.4 


5.7 


3.2 


1.6 



Table 1. Characteristics of the simulated clusters, having T^i^ 3 keV, to 
which we apply the deprojection procedure. Col. 2: virial mass; Col. 3: 
spectroscopic-like temperature; Col. 4; virial radius; Col. 5: rsoo- 



We generate synthetic X-ray and tSZ maps of these clusters, 
to which we apply our method of mass reconstruction. Our proce- 
dure of map making is described in detail in A07, while here we 
summarize the most important characteristics: /) we generate three 
maps of each cluster, by projecting it on the principal axes of the 
inertia tensor, so that the projection on the x, y and z axes represent 
the projection along the minimum, medium and maximum elonga- 
tion of the cluster, respectively; ii) the X-ray maps are obtained 
by Montecarlo generation of 10^ photon counts from the surface 
brightness maps in the [0.5-2] keV energy band, after smoothing 
it with a PSF having a FWHM of 0.5 arcsec (comparable to the 
Chandra one at the aim point), without adding the contribution from 
any instrumental or cosmic background; ;//) the noise setup for the 
tSZ maps is modeled on the Cornell Caltech Atacama Telescope 
(CCAT). The maps are convolved with a Gaussian beam having 
FWHM= 0.44' and then it is added a white noise. In order to have 
approximately the same signal-to-noise for all objects, the noise 
level is 3 /iK/beam for clusters with Tsi < 4 keV and 10 /iK/beam 
for the hotter ones. 

In the following, we will show detailed results on the recovery 
of total mass profiles for a subset of 4 clusters, which are indicated 
in Table [T] with Cl~4. The first three of them are extracted from 
the cosmological hydrodynamical simulation, while C4 is one of 
the massive clusters simulated at higher resolution. C2, C3 and C4 
are typical examples of clusters at low, intermediate and high tem- 
perature, while CI is an interesting case to understand the effect 
of fore-background contaminations. These four clusters have been 
used also in A07 as examples. 



3 THE HYDROSTATIC EQUILIBRIUM 

Total mass profiles from observations of the ICM are computed by 
assuming that the gas lies in hydrostatic equilibrium (HE hereafter) 
within the cluster gravitational potential. In the case of spherical 
symmetry, the equation of hydrostatic equilibrium can be cast in 
the form 

" d In Jie d In T " 



Mtot{< r) = - 



kT{r) r 



dlnr 



dlnr 



(2) 



where T{r) is the temperature at the radius r, /i the mean molecular 
weight (n ~ 0.6 for a gas of primordial composition), and nip the 
proton mass. Note that the mass at a given radius depends only 
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rsoo, a result which is i n line with those found by other analyses of 
simulated clusters (e.g., Nagai et al.l2007l : [Puchwein & Bartelmannl 



0.05 0.1 



r/r . 

/ vir 



Figure 1. The ratio between the hydrcstatic and the true mass profiles, av- 
eraged over the set of simulated clusters. The shaded area encompasses the 
16 and 84 percentiles, while the central continuous line shows the average 
profile of the deviation. 



upon the local pressure derivative and is unaffected by the physical 
properties of the cluster at smaller or larger radii. 

While to first approximation clusters are quite close to the 
condition of pressure equilibrium, small but sizable deviations 
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). These de- 



viations are generally ascribed to the presence of non-negligible 
stochastic gas motions, which provide an effective non-thermal 
pressure support, thereby leading to an underestimat e of the to- 
tal gra vitating mass when not accounted for in eq. iRasia et al] 
| |2004) suggested the addition of an extra-term, which takes into 
account the presence of gas motions, to the equation of the hydro- 
static equilibrium. Jeltema et al. ( 2007) find a correlation between 
the amount of substructures and the underestimate of the total mass, 
in a set of hydrodynamical simulations. However, the large scatter 
in this correlation around the mean relation suggests that substruc- 
tures may not be the only sources of the bias in the mass recon- 
struction. 

Before applying our procedure of mass reconstruction, we as- 
sess the degree of violation of HE in our clusters. This will allow 
us to quantify by how much the differences between reconstructed 
and true mass profiles are due to biases in the deprojection method 
or to violation of the HE. 

We apply eq. ([2} to the true, 3-D density and temperature pro- 
files, by performing a numerical derivative in the log-log space, 
using a 3-point Lagrangian interpolatiorQ. We will refer to the 
mass so computed as the hydrostatic mass, M^yd hereafter. In Fig- 
ure[T]we show the profiles of the ratio Mhyd/ Mtme averaged over 
our set of simulated clusters, along with the Icr scatter around 
the mean. We generally find that the hydrostatic masses under- 
estimate the true ones on average by ~ 10 per cent out to about 



^ The profile of log(nT) is computed for log-equispaced values of r. To 
compute the derivative at a given radius, we fit with the point of interest and 
its two adjacent ones with a quadratic curve. The first-order derivative at the 
middle point is taken to be the derivative of the profile at that radius. 



|2007| ; Ijeltema et al.ll2007l) . At larger radii this underestimate in- 
creases, reaching about 20 per cent at r^ir, along with an increase 
of the scatter. This is consistent with the expectation that outer clus- 
ter regions deviate more from the condition of HE, due to the pres- 
ence of ongoing mergers and continuous gas accretion. 



4 METHODS OF MASS PROFILE RECONSTRUCTION 

The methods described here, aimed at recovering the cluster mass 
profiles, represent a development of the MonteCarlo Markov Chain 
(MCMC) maximum likelihood deprojection technique described in 
A07. The algorithm is modified in order to solve the hydrostatic 
equilibrium equation while deprojecting the cluster images. In this 
way, the temperature profile is computed from the gas density and 
total mass profiles, so that all quantities are derived simultaneously 
and in a fully self-consistent way. We provide here below only a 
brief description of the deprojection technique, while we refer to 
A07 for further details. 

The cluster is assumed to have a onion-skin structure, with 
N concentric spherical shells, having uniform gas properties (gas 
density and temperature). The image of the cluster (in both tSZ and 
X-ray) is divided into A'^ rings, which have the same limiting radii 
of the shells for sake of simplicity. The tSZ and X-ray profiles are 
then computed by projecting this onion-skin model in the plane of 
the sky. We divide the virial radius and rsoo into A^ = 14 and 
A^ = 10 rings, respectively, equally spaced in logarithm. This so- 
lution represents a good compromise between the resolution of the 
profiles and the noise level. Finally, the inner 50 kpc of all clusters 
are excluded by all fits. 

The deprojection is then performed by the maximization of a 
likelihood function, which is computed by comparing the observed 
tSZ and X-ray surface brightness profiles with the ones obtained 
from the model. This approach has the advantage of deprojecting of 
both X-ray and tSZ profiles simultaneously, directly obtaining the 
whole gas density and temperature and total mass profiles, along 
with the corresponding uncertainties. Moreover, it is possible to 
introduce in the likelihood extra terms in order to improve the ac- 
curacy and robustness of the technique. In particular, we adopt a 
regularizatio n constraint, ba sed on the Philips-Towmey regulariza- 
tion method ( iBouchetll 19951 an d references therein). This method 
has been already used also by ICroston et al.l ( |2006|) to deproject 
X-ray imaging and spectral data. It works by interpolating with a 
quadratic curve each group of three consecutive points in the pro- 
file and then minimizing the second derivative (i.e. the curvature) 
of such curves at the position of the middle point (see A07 for a 
more detailed description). The general effect of the regularization 
is to smooth out oscillations in the profiles, which are due either to 
genuine substructures or to noise which propagates from adjacent 
bins in the deprojection. 

We define a joint likelihood for the tSZ profile, Ctsz, and for 
the X-ray surface brightness profile, Cxray, also including a term 
associated to the regularization constraint, -C^eg- Since these three 
terms are independent, the total likelihood is given by the product 
of the individual ones: 



(3) 



The fitting parameters are represented by the gas density and 
total mass profiles (plus the external pressure). The algorithm first 
computes the temperature profile from the density and mass pro- 
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files by inverting the hydrostatic equilibrium equation; then the gas 
density and temperature profiles are combined to compute the pro- 
jected model profiles of X-ray and tSZ; finally these model pro- 
files are compared to those obtained from the mock observations 
to compute the joint likelihood. As for the gas density profile, our 
approach is completely model-independent, since we treat its value 
into each spherical shell as a free parameter, as described in A07. 
As for the reconstruction of the total mass profiles, w e adopt two 
differe nt approaches. Both methods are described in lEttori et al.l 
( l2002h for the application to X-ray datasets, althoug h with a dif- 
ferent implementation. Also. lMorandi & Ettoril jzOOTi) present a re- 
fined version of these methods for t he stu dy of high resolution X- 
ray observations, while lRasia et alj ||2006^ discuss their limitations 
by applying them to a set of hydrodynamical simulations. 

• Method 1. This method consists in the direct numerical inver- 
sion of the HE equation. It does not assume any particular func- 
tional form for the mass profile. Instead, the integrated mass en- 
closed by the mean radius of each shell is treated as a free param- 
eter. The only constraint that we impose is that the mass has to 
increase with radius, in order to avoid unphysical solutions. The 
advantage of this method is that it provides a completely model- 
independent reconstruction of the mass profile, which relies only 
on the assumption of hydrostatic equilibrium. 

• Method 2. This method is based on assuming the functional 
form for the mass profile provided by the NFW model. Its major 
advantage with respect to Method 1 is that the reconstruction be- 
comes more stable, at the cost of assuming a particular model for 
the mass profile. The NFW model is widely adopted in the mass 
reconstruction from X-ray observations. These analyses showed 
that it provides a remarkably good description of the mass pro- 
files of massive ga laxy clusters, out to large portions of their virial 
radii (see also, e.g. Pointe couteau et al.l2005l :IVikhlinin et alj200^ : 
iGastaldello et al...2007.) . The NFW expression for the mass profile 
reads 



M(< r) = 4TvripcritScf{x) 



(4) 



where is a characteristic scale length, x = r /rs is the distance 
from the halo centre in units of r^, f{x) = ln{l + x) — x/{l + x), 
pcrit is the critical cosmic density and 5c is a characteristic over- 
density. It is common (and more convenient) to rewrite the above 
equation by expressing 5c as a function of the concentration pa- 
rameter c = t/s/ts, where A is a given overdensity. Here and in 
the following, we adopt for A the value of the density contrast at 
virialization predicted by the spherical colla pse model, which cor- 
responds to A ~ 100 in our cosmology (e.g.. lEke et aljl996h . The 
characteristic overdensity is then written as 



5c = 



3 f{x = c) ' 

so that the NFW profile becomes 

, ,/ X 47r , , . fix) 

m<r) = -ir.crPcruAjI^^ 



Ma 



(5) 



(6) 



As for the temperature profile, it is then computed in both 
methods from the model mass profile, under the HE assumption 
(see Section[3}. By inverting the HE equation one first obtains the 
gas pressure profile and then the gas temperature, by combining 
pressure with gas density. In this way we obtain 



fcT(r) = 



ne(r) 



G/irrip 



ncM{< r) 



dr - Pijout) \ .(7) 



In the above equation the integral is performed from the outermost 
radius from which the deprojection is performed, rout, and the ra- 
dius of interest r, while Pout is the electron pressure at r = rout- 
In both methods, we apply the regularization constraint to the tem- 
perature profile, in order to smooth out spurious fluctuations which 
are essentially due to the presence of noise. The regularization term 
in the likelihood is given by (see A07 for further details): 



ln{Cr 



JV-l 
i = 3 



(8) 



where A = 2.5 and fi is the value of the temperature in keV in the 
i-th bin. 

Since the hydrostatic equilibrium equation constraints only the 
pressure difference between two poin ts, it is necessary to in troduce 
a further parameter P{rout) (see also lMorandi et al]|2007l) . In par- 
ticular, in the case of Method 1 this parameter is completely degen- 
erate with the mass enclosed in the outermost bin. The mass has 
only a lower boundary, given by the fact that it cannot be lower 
than the mass enclosed by the inner bin. This turns into an insta- 
bility of the fit which generates an overestimate of the global mass. 
This problem is solved by applying the regularization constraint 
to the pressure profile. In this way, in Method 1 the regularization 
constraint is actually applied in two places. Note that once the tem- 
perature profile is regularized, the pressure profile becomes regular 
as well, so its contribution to the global likelihood turns out to be 
generally small. The regularization of the pressure profile is more 
important in the outermost two bins, where it is used to break the 
degeneracy between Pout and the total mass. So, in Method 2 we 
insert in the likelihood a second regularization term, which has the 
same structure of eq. where now A = 5 and fi is the logarithm 
of the pressure in the i-th bin. Note that the numerical value of A 
does not give a direct indication of the weight given to regulariza- 
tion, because this depends on the units used for the quantity fi . 

Given the high number of fitting parameters (2N + 1 for 
Method 1, A'' + 3 for Method 2 respectively), we apply a MCMC 
technique, which is described in A07. The technique has the advan- 
tage of computing the (marginalized) probability distribution of all 
parameters simultaneously, together with the whole covariance ma- 
trix. For a more detailed information on M CM C fitting techniq ues 
we refer to lNeall ( ll993h . lGilks et al.l ( 1 1996b and lMacKavl ( Il996n . 



5 RESULTS 

As already pointed out, both methods of mass reconstruction adopt 
a different implementation with respect to the deprojection tech- 
nique presented in A07. Reassuringly enough, the gas density and 
temperature profiles obtained from all methods are virtually identi- 
cal. For this reason, they are not shown in the following while we 
will concentrate our discussion only on the mass reconstruction. 



5.1 Total mass profiles: method 1 

We show in Figure[2]the comparison between the reconstructed and 
the true mass profiles for the Cl-4 clusters, chosen as examples. We 
generally find an underestimate by about 10-15 per cent through- 
out the virial radius, with slightly larger deviations in the centre 
and in the outskirts. For all these four clusters (and for most of the 
others) the uncertainties in X-ray and tSZ profiles do not allow us 
to place tight constraints on the amount of mass contained within 
the innermost bin. Since this quantity has no lower boundary (e.g. 
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Figure 2. Mass reconstruction for the clusters CI (upper left), C2 (upper right), C3 (lower left) and C4 (lower right), while applying Method 1. The black 
triangles (squares, circles) and line represent the reconstruction along the x {y, z) axis. The red dashed line represents the true mass profile, while the blue 
dotted line is hydrostatic mass obtained from the application of eq. (2) to the true 3-dimensional gas density and temperature profiles. The lower part of each 
panel shows the ratio between the same quantities and the true mass. 



given by the mass contained within an inner bin), we find cases in 
which the lower limit of this mass takes very low or even nega- 
tive values. To avoid this problem, we fix by hand a lower limit for 
the mass within the innermost bin at 10^^ h~^MQ. A comparison 
of the reconstructed mass profiles with the profiles of hydrostatic 
mass, Mhyd, generally show a close agreement. This demonstrates 
that the main source of systematics is intrinsic, namely the devi- 
ation of the gas from a perfect hydrostatic equilibrium, while the 
deprojection method is essentially unbiased. While this is strictly 
true out to r ~ 0.5r„ir, we note a tendency for the reconstructed 
mass to lie above Mhyd when approaching the virial radius. We 
attribute this to the lower signal-to-noise ratio in the external re- 
gions (SNR ~ 4 — 5 for both the SZ and X-ray signals in the 
outermost bin, compared to 15-20 in the innermost bins) and to a 
larger impact of fore/background contaminations in these regions. 
In fact, the recovered mass at these radii is bound not to be smaller 
than the mass at the inner radii, while it does not have any upper 
boundary. As a consequence, any deviation associated to noise or 
contamination can only act in the direction of increasing the mass. 



We further point out that our synthetic X-ray maps do not include a 
background, so the noise here is only given by the poissonian noise 
of the signal. Clearly, the inclusion of a realistic background would 
limit our analysis to scales comparable to those reached by X-ray 
observa tions, which only in some cases trace the emissivity beyond 
~ r2oo ( lNeumanj|2005l) . 

These results are also confirmed from the analysis of the 
whole set of simulated clusters. Indeed, Figure[3]shows that in the 
inner regions (out to about r = 0.15r^ir) the mass profile recov- 
ered from our deprojection algorithm are quite close on average 
to the hydrostatic masses, while being consistently smaller that the 
true mass profiles, by about 10-15 per cent. Instead, in the outskirts 
the reconstructed mass profiles seem to well recover the true ones. 
However, this agreement is due to the effect of two biases, namely 
a violation of HE and the effect of noise in the deprojection algo- 
rithm, which act in opposite direction, rather than to an intrinsically 
better performance of the deprojection method at large radii. 

In A07 we showed that the CI cluster has a significant contam- 
ination in the projection along the z axis, due to a merging group 



Mass profiles from SZ/X-ray in simulated clusters 



1.5 - 



0.5 




0.5 - 



0.05 0.1 



0.5 



r/r . 

/ vir 



Figure 3. Reconstructed vs. trae (upper panel) and hydrostatic (lower panel) 
mass, averaged over the set of simulated clusters, while applying Method 
1. In each panel, the horizontally (vertically) shaded area represents the 
meanitlcr over the projections along the z (x and y) axis. The blaclc line 
represents the mean over all the projections of all clusters. 



which is approaching the cluster virial region. The effect of this 
contamination was to boost the temperature profile by about 20 per 
cent. In the upper left panel of Figure |2] we show that this causes 
a mass overestimate by a comparable amount of the mass recon- 
structed from the projection along the z axis, which lies above the 
hydrostatic mass in this particular case. Also for the C2, C3 and C4 
clusters the mass reconstructed from the projection along the axis 
of maximum elongation is larger than the mass recovered from the 
other two projections. Figure[3]confirms that this behaviour is still 
present when averaging over the whole set of simulated clusters. 

5.2 Total mass profiles: method 2 

We show in Figure|4]the mass profiles of the same 4 clusters chosen 
as examples, as reconstructed from the application of the method 
based on assuming an NFW model. In this case, the reconstructed 
profiles are not reported with errorbars, which are assigned to the 
total mass M and to the concentration parameter c, which are the 
only two parameters entering in the fitting procedure. The accuracy 
in recovering these parameters is discussed in the following. 

Clearly, assuming an analytical form for the mass profiles has 
the advantage of providing a much more stable reconstruction and 
a smooth profile. Similarly to the case of the Method 1, the recov- 
ered mass is much closer to the hydrostatic mass than to the true 
one. The left panel of Figure [5] shows the accuracy of the recov- 
ered mass profile, after averaging over the whole set of simulated 
clusters. From the upper panel of this plot, we note that the mass 
profile is generally underestimated by about 15-20 per cent. The 
slightly larger underestimate with respect to Method 1 is due to 
the fact that here we are adopting an analytical model. We remind 
that the violation of hydrostatic equilibrium is larger at r„i,., which 
is the outermost radius from which the reconstruction of the mass 
profile begins. Since we force the profile to follow the NFW shape 
when reconstructing at smaller radii, this larger mass underestimate 
is now propagated inward. This is different from what happens for 
the Method 1, where the value of the mass at an inner radius is not 
forced by the extrapolation of the profile recovered at a larger ra- 



true 



'"A. /it 

A.,true 



MA.true 



0.81 ± 0.10 



0.91 ±0.05 0.77 ±0.11 



r < r-500 500 
vir 



0.89 ± 0.16 
0.90 ±0.21 



0.95 ±0.08 
0.95 ±0.10 



0.87 ±0.20 
0.89 ±0.27 



Table 2. Summary of the results of Method 2 over the whole set of simu- 
lated clusters. The first column reports the limiting radius of the fit, while 
the second column reports the overdensity A at which all quantities are 
computed. Columns 3 to 5 report the mean and standard deviation for the 
accuracy in recovering the following quantities: the total mass within Ra 
when it is computed from simulation data; the value of Ra obtained from 
the fit; the total mass enclosed by Ra when it is fitted from data (note that 
in this case the fitted mass and the true mass refer to different radii). All 
these data are computed by excluding the cluster labeled C5 (see text for 
details). 



dius. The right panel of Figure|5]shows the accuracy in recovering 
the total mass as a function of the true mass of the cluster. We find 
that the statistical uncertainties in the estimate of the cluster masses 
are generally smaller than the bias induced by the violation of the 
hydrostatic equilibrium and are also smaller than the scatter associ- 
ated to the choice of the projection direction. This panel also shows 
the presence of an outlier, which is the cluster labeled C5; this ob- 
ject has a much more irregular structure when compared with the 
rest of the sample. For this reason, here and in the following we 
exclude it from the computation of the mean and the standard devi- 
ation over the sample of simulated clusters, as reported in Table |2] 
Within Tjjir, we find Mrec/Mtrue = 0.81 ± 0.10, with no signifi- 
cant dependence on the cluster mass. 

The increase of the HE violation at large radii suggests that 
the mass profiles should be better recovered in normalization by 
limiting the fit to a smaller region, typically ^ r^oo- This radius 
represents the outermost limit typically reached by X-ray observa- 
tions. We repeat the analysis by limiting the fit to rsoo, while using 
the recovered NFW model to extrapolate the mass profile out to the 
virial radius. The left panel of Figure|6]demonstrates that restricting 
the reconstruction of the profile within rsoo has the expected effect 
of reducing the mass underestimate. While the recovered mass pro- 
files are very close to the hydrostatic ones out to the fitting radius 
rgoo, it exceeds the HE predictions at larger radii, while remaining 
within the 10 per cent underestimate of the true mass out to the 
virial radius. This demonstrates that the NFW model reconstructed 
within rsoo can be safely extrapolated to larger radii to extend the 
mass profile reconstruction. After averaging over the set of simu- 
lated clusters, we find that M500 and M„ir are both recovered with 
an underestimate of only 10 per cent (see Table|2j. The right panel 
of Figure[6]shows the accuracy in recovering the virial mass of each 
cluster. This plot has a larger scatter with respect to Figure [5] due 
to the fact the we are now fitting over mass profiles over a nar- 
rower radial range (see Table [2). In agreement with the results of 
IPiffaretti & Valdaminil ilOOA we conclude that a reliable proce- 
dure is to reconstruct the mass profiles within relatively small radii, 
~ '"500, where the HE is not seriously violated, while extrapolating 
to so-obtained mass profile according to the NFW model. 

In the results shown so far, we assumed that we exactly know 
the value of the radius, either r^ir or rsoo, within which the masses 
are recovered. In the analysis of observational data, instead, the 
virial radius is generally not known in advance, but it is estimated 
directly from the recovered mass profile. Therefore, if the mass pro- 
file is under/overestimated, also the es timate of the virial radius will 
be biased low/high. As pointed out bv lNagaietal] ( l2007t) . this turns 
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Figure 4. The reconstruction of the mass profiles for the CI, C2, C3 and C4 clusters, based on the Method 2, compared with the true profiles (dashed curves) 
and with the profiles recovered from the hydrostatic equilibrium (dotted curves). For each cluster, the three solid curves show the reconstruction along the 
three axes of projection. 



out into a larger bias in the recovered Mr,oo or M^ir- In order to 
quantify this effect, we show in the left panel of Figure|7]the virial 
radii inferred from the reconstructed mass profiles, Vvirju, vs. the 
true ones rvir,3im - Clearly, the recovered virial radii are generally 
smaller than the true one, by about 10 per cent. This bias decreases 
by a few percent if r-oir.fu is recovered from the extrapolation of 
the profile fitted within rsoo (seeTable|2l(. The right panel of Fig- 
ure |7] shows the recovered virial mass when using r^irj'u versus 
the true virial mass. A comparison of this plot with the right panel 
of Figure [5] shows that the effect of using the recovered virial ra- 
dius is that of slightly increasing the underestimate of the virial 
mass, which is now Mvir,fit/Mvir,true = 0.76 ± 0.11. Again, 
we verified that this effect is reduced, with M^irjit/Mvir,true = 
0.89 ± 0.27 when using instead the virial radius extrapolated from 
the analysis performed within rsoo . 

Finally, we find that observing the same cluster along different 
lines of sight cause differences in the reconstructed mass by only 
a few per cent (see Figures |5] [6] and [7}. It is worth noting that the 
axes of projection are not chosen randomly, instead they are fixed 
to the principal axes of the inertia tensor. This choice maximizes 



the difference between the projected images, so our estimate can 
be considered an upper limit of this effect. 



6 THE CONCENTRATION-MASS RELATION 

In this section, we first study the c~M relation on the three- 
dimensional mass profiles computed directly from the simulation 
data, so that no observational effects are included. Then, we will 
show how well this relation can be reconstructed from the analy- 
sis of combined X~ray/SZ observations (using Method 2), for the 
subset of clusters having Tsi<^ 3 keV. Since the behaviour of the 
c — M relation strongly depends on cosmological parameters, we 
limit our analysis to those objects belonging to the first simulation 
(i.e. CI, C2, C3 and from C5 to CI 1). Finally, we also com pare our 
results with the mass-concentration relation proposed by E ke et al.| 
( l200lh . using the cosmological parameters assumed in our simula- 
tions (see Section 2). 

Figure[8]shows the c-M relation computed over the 1 17 clus- 
ters identified in the simulation box and having mass assigned by 
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Figure 5. Left panel: the ratio of the reconstructed mass profiles, based on Method 2 (fitting out to r^ir), with the true (upper panel) and hydrostatic (lower 
panel) mass profiles, after averaging over the set of simulated clusters. The shaded area encompasses 68 per cent of the recovered profiles, while the solid line 
shows the median profile. Right panel: the ratio between the recovered and the true masses, both computed within r^ir, as a function of the true mass of each 
cluster. Diamonds, triangles and stars are for the reconsfi'uction from the projection along the x, y and z axis respectively. Errorbars represent the Ic statistical 
uncertainty on the recovered mass, due to the noise in mock SZ/X-ray images. The horizontal dotted line indicates the mean deviation of the reconstructed 
mass. 
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Figure 6. The same as in Figure[5] but fitting the NFW profile out to rsoo, and extrapolating it out to r^i^. 



a frie nds-of-friends algorithm above 10^'' h ^Mq jBorgani et al.l 
l2004h . The left panel of the figure shows the result of fitting the 
integrated mass profile to the NFW formula out to r^jir, while in 
the right panel we exclude the central region with r < O.OSr^i,-- 
Note that excluding the inner regions from the fit has the effect 
of decreasing the resulting va lue of c, thus brin ging the relation in 
closer agreement with that bv lEke et al.l ( 1200 1[) . which is based on 
purely N-body simulations. Including gas cooling and star forma- 
tion leads to the formation of a stellar component, which generally 
has a more sharply peaked mass distribution with respect to the 
dark matter one. As a consequence, one observes a steepening of 
the total mass profile in the core, with a subsequent increase of the 
concentration parameter. Furthermore, gas cooling is also known 
to induce adiabatic contraction, i.e. a steepening of the dark matter 



profile in rea ction to the strong i ncrease of gas density in the core 
regions (e.g. iGnedin et"ID l2004l). We point out that the increase 
of c is much reduced when fitting the differential mass profiles, 
which gives less weight to the inner halo regions, rather than the 
integrated mass profile. However, the case of fitting the integrated 
mass is more relevant for tSZ and X-ray studies, since the appli- 
cation of the hydrostatic equilibrium equation gives the integrated 
mass profile, rather than the differential one. 

In Figure [S] the colour of each symbol indicates the 
spectroscopic~like temperature of the corresponding cluster. We di- 
vide the clusters into 3 temperature bins. We find that hotter clus- 
ters tend to have a slightly larger c than colder ones, at a fixed mass. 
The effect can be seen more clearly in Figure |9l which reports the 
fractional deviation of the spectroscopic-like temperature, ATsi, 
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Figure 7. Left panel: the estimated virial radius r^ij-jn using Method 2 out to r^ir) vs. the true one computed from the simulation data r^ir^sim- The dotted 
line reports the one-to-one relation. Right panel: the recovered mass within r^irju vs. the true mass within r^ir^sim- In both panels, diamonds, triangles 
and stars are for the reconstruction from the projection along the x, y and z axis respectively. Errorbars represent the Icr statistical uncertainty on the recovered 
mass, due to the noise in mock SZ/X-ray images. The horizontal dotted line represents the mean over all the mass determinations. 




Figure 8. The concentration-mass relation for all the clusters extracted from cosmological box o f lBorganiet^j2004) . Diamonds, triangles and circles are for 
objects with Tgi < 1.5 keV, 1.5 < Tgi < 2 keV and Tgi > 2 keV, respectively. The objects used i n this paper are i ndicated by filled symbols. In both panels 
the dotted line represents the least-squares fit on the points, the solid line is the model proposed bv lEke et alj j200il) . The left and the right panels correspond 
to fitting over the whole radial range, out to r„i,., and after excluding the central region within 0.05r„i^, respectively. 



and concentration Ac from the mean Tsi-M and c-M relations 
respectively. We find that the correlation factor is quite high, with 
r = 0.86. This correlation is explained by considering that clus- 
ters with larger c have a deeper gravitational potential and thus a 
hotter ICM. We point out that this does not depend on the use of 
the spectroscopic-like temperature. In fact, when using the mass- 
weighted temperature, Tmw = X/ ^iTi/ mi, we find that this 
correlation is only slightly weaker, with r = 0.77, but still signifi- 
cant. As a consequence, selecting the clusters through their (X-ray) 
temperature (as we do in this paper) has the effect of biasing high 
the recovered c-M relation. 

Our resul t s go i n the opposite direction with respect to those 
by I Yang etal] j2008l) . who instead find a negative correlation be- 



tween temperature and concentration in a set of simulated galaxy 
clusters. These authors based their analysis on a cosmological box 
simulated with the Eulerian FLASH code. Tracing the origin of this 
difference between our and their results would require a detailed 
comparison. We note that the run analysed by Yang et al. 1 1 I2OO8I) 
includes quite difference physical processes: it assumes gas pre- 
heating at z = 3, while it does not include radiative cooling, star 
formation and stellar feedback. We also note that the size of the 
grid in their simulation is 250 /i^^kpc, thus implying that rsoo in a 
typical cluster is sampled only with few resolution elements. 

The results discussed in this paragraph and shown in Fig- 
ures [TO] and [TT] have been obtained by applying Method 2 out to 
r^ir- The procedure is discussed in detail in Section [J!2l In Fig- 
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Figure 9. The correlation between the scatter in the temperature-mass re- 
lation and the scatter in the concentration-mass relation. Results are shown 
for the same objects as in Fig. [8] 
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Figure 11. The difference between the recovered and the true concentration 
parameter as a function of the true cluster mass within r^ir- The errorbars 
show the 1 (T statistical error from the recovery of the cluster mass profile. 
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Figure 10. The concentration-mass rel ation recovered from deprojection. 
The solid line represents the model by Eke et alj jlQOlh . The ellipses show 
the la confidence regions. Diamonds, triangles and stars are for the recon- 
struction from the projection along the x, y and z axis respectively 



ure[TO]we show the c-M relation, as it would be recovered from the 
deprojection of SZ/X-ray images. We find that it is systematically 
higher than the expected one. This shift is due to the temperature 
cut in selecting clusters, which tends to favor more concentrated 
systems. Finally, we show in more detail in Figure[TT]our results on 
the concentration parameter. We report the error on c as a function 
of the cluster mass. We find that (crec — ctme) ~ 0.1 ± 1.5, with 
no significant trend with the cluster mass. This result demonstrates 
that the value of c is recovered unbiased on average, although with 
a fairly large scatter. 



Correctly measuring the total collapsed mass of galaxy clusters is 
of fundamental importance for these objects to be used as tools for 
precision cosmology. In cluster studies based on the observations 
of the intra-cluster medium (ICM) the mass is inferred by applying 
the equation of hydrostatic equilibrium (HE) under the assumption 
of spherical symmetry, thus requiring that both the gas density and 
the temperature profiles are measured to good accuracy. The ob- 
servational analyses carried out to so far are essentially based on 
X-ray data. Although X-ray surface imaging provides robust mea- 
surements of the gas density profiles out to a significant fraction 
of the cluster virialized regions, temperature measurements from 
spatially resolved X-ray spectroscopy are in general much noisier 
and restricted to smaller radii. In this paper we addressed the issue 
of measuring cluster masses by using a combination of Sunyaev- 
Zeldovich (SZ) and X-ray imaging, thereby avoiding X-ray spec- 
troscopy. 

To this purpose we presented in this paper a development 
of the maximum-like lihood deprojection technique described by 
lAmeglio et ai] ( |2007|) , in which we implemented the solution of 
the hydrostatic equilibrium equation, so as to derive profiles of 
gas density, temperature and total mass simultaneously. We applied 
this method to cosmological hydrodynamical simulations of galaxy 
clusters. After quantifying the intrinsic level of HE violation in our 
simulations, we applied two different methods to recover the total 
mass profile from the deprojection of SZ and X-ray maps, both 
methods being based on the assumption of HE. Method 1 adopts 
a model-independent approach, in which the fitting parameters are 
represented by the values of the gas density and total mass profiles 
within 15 radial bins, with the only constraint of increasing mass 
with r adius. Method 2 is in stead based on assuming the mass pro- 
file by iNavarro et ai] ( 1 19971) (NEW). In this case, the mass profile 
is characterized by only 2 fitting parameters, namely the concentra- 
tion c and the total mass M. 

The main results of our analysis be summarized as follows. 

Rasia et al.l 



In keeping with previous analyse s (e.g. 



2004 lKavetal.1 |2004 iRasia et al.l l2006l : iNagai et alj |2007|: 
Piffaretti & Valdaminil 20081) , we find that deviations from hydro- 
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static equilibrium are quite common in simulated clusters, due 
to the presence of pressure support associated to stochastic and 
bulk gas motions. These deviations lead to an average mass un- 
derestimate of about 10 per cent within rsoo- It increases at larger 
radii, reaching about 20 per cent at r-uir, due to the presence of 
ongoing mergers and continuous gas accretion, which characterize 
the cluster outskirts. 

• From Method 1, we find that in the inner regions (out to about 
r = 0.15ri,ir) the mass is recovered with an underestimate of about 
10-15 per cent, largely due to the violation of hydrostatic equilib- 
rium. On the other hand, the reconstructed mass becomes closer 
to the true one when moving to the outskirts. In this case, the ef- 
fect of noise in the deprojection algorithm tends to compensate the 
underestimate due to the violation of the hydrostatic equilibrium. 
Generally, the fairly large number of fitting parameters causes a 
significant scatter (about 15 per cent) in the recovered mass pro- 
files. 

• As for the Method 2, we find that the total mass is underesti- 
mated by a larger amount, about 20 per cent, when the deprojection 
is performed from r^ir- Indeed, using a fixed NFW functional form 
for the mass profile forces the mass underestimate found at the out- 
ermost radius to propagate to smaller radii. Therefore, the overall 
lower normalization of the mass profile is just the consequence of 
the larger violation of the HE found at r^ir- On the other hand, as- 
suming a fixed functional form for the mass profile leads to a more 
stable, although more biased, reconstruction, with a scatter of about 
10 per cent, thus lower than found with Method 1. 

• As expected, fitting the model profiles in the deprojection 
within rsoo and extrapolating out to the virial radius is a safer pro- 
cedure, thus in agreement with the results by Piffaretti & Valdarnini 
( |2008[) . The mass underestimate is reduced to ~ 10 percent, but at 
the expense of increasing the cluster-to-cluster scatter to about 20 
per cent. 

• We verified that the choice of the line of sight generally affects 
the mass reconstruction. With Method 1, we show that the mass 
reconstructed from the projection performed along the maximum 
elongation axis is generally larger, by about 10 per cent, than from 
the other two projection directions. Using instead Method 2, we 
find that the typical scatter between the masses reconstructed from 
different projections of the same cluster is typically of only a few 
percent. 

• The relation between concentration parameter of the NFW 
density profile and cluster mass shows a strong correlation with 
the temperature of the cluster: at a fixed mass, hotter clusters tend 
to be more concentrate. For this reason, in our selection of clus- 
ters having T^i ^ 3 keV we find a c — M relation having higher 
norm alization than that calibrated from N-body simulations (e.g., 
lEke et al. 2001). 

• Using Method 2, the concentration parameter of the NFW pro- 
file is recovered on average without any significant bias, but with a 
significant scatter of Ac ~ 1.5. 

Our results lend support to the efficiency of combining X-ray 
and tSZ imaging data to recover the total mass profiles of galaxy 
clusters. In fact, the main bias that we found is intrinsic, since it 
is due to deviations from perfect hydrostatic equilibrium. This ap- 
proach has several advantages with respect to the traditional one 
based on X-ray spectroscopy. Firstly, the temperature recovered 
from the fit of the X-ray spectra generally differs from the electron 
temperature, by an amount which depends on the degree of com- 
plexity of the I CM thermal structure (e.g., iMazzotta et alj|2004l : 
IVikhlininlliooeb . Secondly, X-ray surface brightness profiles can 



be obtained with good precision with a relatively small number of 
photon counts (~ lO'^), while at least ten times more photons are 
required for a reliable spatially resolved spectroscopy. Also, once 
the cosmic and instrumental backgrounds are under control, this 
opens the possibility of tracing the surface brightness over a large 
portion of the cluster virial regio ns, as already demons t rated with 
ROSAT-PSPC imaging data (e.g.. lVikhlinin et alJll999l : lNeumannl 
2005). Since the tSZ has the potential of covering a large range in 
gas density, its combination with low-background X-ray imaging 
data would allow one to better characterize the outskirts of galaxy 
clusters. Moreover, another advantage is that the tSZ is indepen- 
dent of redshift, so its combination with X-ray surface brightness 
will substantially improve the analysis of cluster properties at high 
redshift, where accurate X-ray spectroscopy may be very hard, if 
not impossible, to obtain. 

A limitation of the analysis presented in this paper is that we 
did not include realistic backgrounds in the generation of the X- 
ray and tSZ maps. As we have just mentioned, there are interest- 
ing perspectives for a good characterization of the X-ray back- 
ground. However, the situation could be more complicated for the 
tSZ background. In this case, contaminating signals from unre- 
solved point-like radio sources (e.g., Bartlett & Melin 2006) and 
fore/background galaxy groups (e.g.. iHallman et al.l 120071) could 
affect the tSZ signal in the cluster outskirts. In this respect, the pos- 
sibility of performing multi-frequency observations with the good 
angular resolution offered by interferometric techniques will surely 
help in characterizing and removing these contaminations. 

Single-dish sub-millimetric telescopes of the next generation 
promise to provide tSZ images of clusters with a spatial resolution 
of few tens of arcsec, while covering fairly large field of views, 
with 10-20 arcmin aside, with a good sensitivity. At the same time, 
future satellites for X-ray surveys (e.g. eROSITA) will have the ca- 
pability of surveying large areas of the sky with high sensitivity and 
good control of the background. These observational facilities will 
open the possibility of carrying out in survey mode high-quality 
tSZ and X-ray imaging for a large number of clusters. The appli- 
cation of deprojection methods, like those presented in this paper 
will provide reliable determinations of both the gas mass and total 
mass profiles. This will greatly help to fully exploit the potentiality 
of galaxy clusters as tools for precision cosmology. 
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